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Abstract 



The nature of the fuel that drives today's cosmic acceleration is an open and 
tantalizing mystery. We entertain the suggestion that the acceleration is not 
the manifestation of yet another new ingredient in the cosmic gas tank, but 
rather a signal of our first real lack of understanding of gravitational physics. 
By requiring that the underlying gravity theory respects Birkhoff's law, we 
derive the modified gravitational force-law necessary to generate any given 
cosmology, without reference to the fundamental theory, revealing modifi- 
cations of gravity at scales typically much smaller than today's horizon. We 
discuss how through these modifications, the growth of density perturbations, 
the late-time integrated Sachs-Wolfe effect, and even solar-system measure- 
ments may be sensitive to whether today's cosmic acceleration is generated by 
dark energy or modified gravitational dynamics, and are subject to imminent 
observational discrimination. We argue how these conclusions can be more 
generic, and probably not dependent on the validity of Birkhoff's law. 
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I. INTRODUCTION 



The discovery of a contemporary cosmic acceleration [1,2] is one of the most profound 
scientific observations of the latter part of the last century. What drives that acceleration 
remains an open and tantalizing question. A "conventional" explanation exists for the cause 
of that acceleration: vacuum energy provides the necessary repulsive gravity in general 
relativity to drive an accelerated expansion of the universe. Variations on this vacuum- 
energy theme, e.g. quintessence, promote the energy density to the potential energy density 
of a dynamical field. Such additions to the roster of cosmic sources of energy-momentum are 
collectively referred to as dark energy. Accounted for as such, dark energy would constitute 
the majority of the energy density of the universe today. 

However, what if one were to view the current cosmic expansion not as yet another in- 
gredient in our already complex cosmic soup, but rather as a signal of our first real lack 
of understanding of gravitational physics? (Although others have argued that dark matter 
was that first signal of post-Einsteinian gravity [3,4], we will not address ourselves here to 
that possibility.) An instructive example along that direction is the braneworld theory [5] of 
Dvali, Gabadadze, and Porrati (DGP). In this theory, gravity appears four-dimensional at 
short distances but is altered at large distances through the slow evaporation of the graviton 
off our four-dimensional braneworld universe into an unseen, yet large, fifth dimension [5-7] . 
DGP gravity provides an alternative explanation for today's cosmic acceleration [8,9]: just as 
gravity is conventional four-dimensional gravity at short scales and appears five-dimensional 
at large distance scales, so too the Hubble scale, H{t), evolves by the conventional Fried- 
mann equation at high Hubble scales but saturates at a fixed value as H(t) approaches a 
value equivalent to the inverse of the crossover distance between four and five-dimensional 
behavior, r^^. Thus, if one were to set that crossover distance scale to be on the order of 
Hq^, where Hq is today's Hubble scale, DGP could account for today's cosmic acceleration 
in terms of the existence of extra dimensions and a modifications of the laws of gravity. 

We would naively expect not to be able to probe this extra dimension at distances much 
smaller than the crossover scale Vc = Hq^. However, in DGP, although gravity is four- 
dimensional at distances shorter than Tc, it is not four-dimensional Einstein gravity - it is 
augmented by the presence of an ultra-light gravitational scalar. One only recovers Einstein 
gravity in a subtle fashion [10-13], and a marked departure from Einstein gravity persists 
down to distances much shorter than Vc- For example, for Tc = Hq^ and a central mass 
source of Schwarzschild radius Vg, significant and cosmologically-sensitive deviations from 
Einstein gravity occur at distances greater than [12-15] 




(1.1) 
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Thus a marked departure from conventional physics persists down to scales much smaller 
than the distance at which the extra dimension is naively hidden, or for our discussion 
here, the distance at which the Priedmann equation was modified to account for accelerated 
cosmic expansion. (Other theories have since shown how alternative modifications of gravity 
at large distances can lead to late-time acceleration without dark energy [16,17].) 

Recently, appeals have been made to a direct empirical modification of the Priedmann 
equation in order to explain cosmic acceleration without dark energy [18-20], broadening 
the theme explicitly and completely realized by DGP braneworlds. Unfortunately, fully 
self-consistent models implementing these ideas have yet to be found that can reproduce 
the desired general modifications to the Priedmann equation, making it difficult to establish 
the full consequence of the proposed new physics. An important question to be asked is 
whether gravity theories that produce such modified Priedmann equations can be devoid of 
observable consequences other than late-time acceleration? There may be something to be 
learned from the example of DGP gravity. Will these new theories generically lead to similar 
deviations from Einstein gravity at scales much smaller than today's Hubble radius? 

Pursuing this line of thought without an underlying model is difficult, but not altogether 
impossible. In this paper, we allow ourselves an assumption about the structure of a possible 
modified theory of gravity, taking the constraint of Birkhoff 's law as the example, and show 
how with that assumption alone, one may kinematically ascertain gravitational force inter- 
actions, avoiding reference to fundamental dynamics. We then describe how modification of 
cosmology at today's Hubble scale can naturally affect gravitational physics at much smaller 
(e.g. astrophysical, and even under certain circumstances solar system) distance scales. It 
is worth noting that DGP gravity does not respect Birkhoff's law, and we conclude with a 
discussion on how short-distance modifications of gravitational interactions may be a general 
consequence of gravity modified at cosmological scales. 

II. THE FRIEDMANN EQUATION AND THE FORCE LAW 

The premise of our approach is that the static, Schwarzschild-like metric, or more specifi- 
cally its geodesies, may completely determine cosmological evolution, i.e., that the cosmolog- 
ical evolution is driven not by some dark energy background, but by the matter content itself. 
This should be a familiar premise - it underlies the usual undergraduate-level derivation of 
the Priedmann equation from Newton's laws of motion and his Law of Gravitation. One 
posits that in a homogeneous isotropic dust-filled universe, the acceleration of a comoving 
(geodesic) observer a distance r from some arbitrary origin, is determined by the gravita- 
tional pull of the matter interior to the sphere of radius r. One then derives an equation for 
the evolution of r with time. That one arrives at the precisely correct General-Relativistic 
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Priedmann equation, may seem surprisingly coincidental but is really a consequence of homo- 
geneity, energy conservation and Birkhoff 's law. Thus in the standard cosmology, Newton's 
Law of Gravitation, or now its General Relativistic generalization, and the observed scale 
factor evolution are used to extract the stress-energy content as a function of time. The 
approach here will be to fix the form of the stress energy to be that of dust (non-relativistic 
matter) and use the observed evolution of the scale factor to extract a new generalization 
of Newton's Law. 

We take the cosmological metric to be the Robertson- Walker metric (flat space, for 
convenience) 

ds^ = - a^{t) dijdx'dx^ 

= dt^ - a^{t) [dX" + XMq\ . (2.1) 

We imagine then that we are given a complete cosmological evolution, i.e., a specified cos- 
mological scale factor evolution, a{t). (The understanding being that a{t) is determined on 
the past light cone by observations, and that the assumption of spatial homogeneity carries 
this to the whole space.) If we were to solve the Einstein equations for this metric, we would 
be forced to specify a matter content and an equation of state as a function of time. Such 
an equation of state would not normally be matter dominated, and we would be forced to 
consider the non-dust component to be a sort of dark energy. However, we will instead 
presume that the scale factor a{t) is the result of a pure dust configuration, and ask what 
modifications of gravity would be necessary to yield the given a{t). 

Such a prescription for modifying the Einstein equations is not unique. However, we can 
require that this hypothetical new gravity theory obeys a generalization of Birkhoff 's law: 
for any test particle outside a spherically symmetric matter source, the metric observed by 
that test particle is equivalent to that of a point source of the same mass located at the 
center of the sphere. With that one restriction on the theory, we can deduce completely the 
Schwarzschild-like metric of the new hypothetical gravity theory that gives us the prescribed 
cosmological evolution with a dust-filled universe. 

The procedure for determining the metric from the cosmological evolution is as follows. 
Consider a uniform sphere of dust. Imagine that the evolution inside the sphere is exactly 
cosmological, while outside the sphere is empty space, whose metric (given Birkhoff 's law) is 
Schwarzschild-hke (as defined by the metric Eq. (2.2) below). The mass of the matter source 
(as determined by the form of the metric at short distances) is unchanged throughout its 
time-evolution. The surface of the spherical mass therefore charts out the metric through 
all of space as the sphere expands with time, so long as we demand that the cosmological 
metric just inside the surface of the sphere smoothly matches the Schwarzschild solution 
just outside. In order to see how the metric depends on the mass of the central source, we 
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just take a sphere of dust of a different initial size, and watch its surface chart out a new 
metric. 

We start with a Schwarzschild-hke metric in the usual form 

^ goo{r)dT^ - grr{r)dr^ - r^dn , (2.2) 

and rewrite it in the new form 

ds^ = N\t, \)de - a^{t)dX^ - r^dfl , (2.3) 

where now r = r{t, A) and a{t) is the given cosmological evolution. In order to determine 
the coordinate transformation T = f{t,X) and r(t, A), we equate the forms Eq. (2.2) and 
Eq. (2.3): 

^ goof' - 9rrr' (2.4) 
= gooff - grr-rr' (2.5) 
-a' = goof" - grrr" , (2.6) 

where dot denotes partial differentiation with respect to t (holding A fixed) while prime 
denotes partial differentiation with respect to A (holding t fixed). 

The object of transforming to the form Eq. (2.3) is that the bounding surface of the 
spherical mass distribution can be taken to be at fixed A = A*. We wish now to make an 
identification between the interior metric Eq. (2.1) and the exterior metric Eq. (2.3) which is 
smooth, and such that this matching surface is also a geodesic of the exterior (Schwarzschild- 
like) metric. This is realized if the following conditions hold at the boundary: 



r{t,X^) ^ X^a{t) (2.7) 

r'(t. A.) = a(t) , (2.8) 

N{t,X,)^l (2.9) 

7V'(t,A,)=0. (2.10) 



Using Eqs. (2.7) and (2.8) as the boundary condition, one can integrate Eqs. (2.5) and (2.6) 
to arrive at the complete coordinate transformation T — f{t,X) and r{t,X), for arbitrary 
functions g'oo('^) and grr{r)- For Eqs. (2.9) and (2.10) to be satisfied, conditions need to be 
placed on g'oo(^) and grr{f)'- 

r\t,X,) = l-g-\r{t,X,)) (2.11) 
googrr — E'^ — constant , (2-12) 
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for Eqs. (2.9) and (2.10) respectively. One can quickly confirm that these expressions com- 
bined imply A = A* follows a geodesic worldline and, moreover, these expressions allow us 
to determine the metric components of Eq. (2.2) uniquely from a given a{t): 

goo = E\l - Xy) (2.13) 
^-^-l-A^a^, (2.14) 

with r(t, \^) = A*a(t). 

Thus, by requiring our new gravitational physics to respect BirkhofF's law, the metric 
around a spherically symmetric matter source is completely specified by the cosmology 
a{t). To determine the dependence of the metric on the source mass, one need only select 
a different A*. The remaining parameter E may be specified by an arbitrary choice of 
time normalization. Notice that for a general scale factor evolution a{t), superposition and 
linearity of the metric must be sacrificed, even in the weak-field limit. ^ Let us see what 
consequences this has for phenomenology. 



III. GOVERNING SCALES 



The exterior (Schwarzschild-hke) metric is given by Eqs. (2.13) and (2.14). For the 
gravitational force law to approach Einstein's at short distances, we require 

Wr)=^-/(r) = l-^, (3.1) 

when r is small and Vg = 2GM is the usual definition of the Schwarzschild radius. On the 
other hand, cosmology at early times (but still during matter dominated regime) must obey 
the conventional Friedmann equation 

/f^-^ = ^P. (3.2) 
3 

where we use p{t) to denote the matter density. The relationship between mass and A* is 
therefore 

M{X.) ^ ^ = ^p{ty = ^XW , (3.3) 



Note that DGP gravity docs not satisfy the condition that gooQrr = constant and, therefore, 
cannot satisfy a dynamical version of Birkhoff 's law. But, while DGP gravity does not fall under 
the category of modified-gravity theories under consideration, it does share many of the same 
properties, such as the lack of superposition and linearity, as well as properties to be noted in the 
coming sections. 
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where p{t) is the matter density, and pa^ is a constant with respect to t for a dust-filled 
universe. 

How small should r be, or how large should H be, for this conventional behavior 
(Eqs. (3.1) and (3.2)) to be apphcable? Clearly, the cosmological evolution must be conven- 
tional when H ^ Hq, where Hq is today's Hubble scale (approximately the scale at which 
acceleration sets in) or in other words, rearranging Eq. (2.14), when 



1 - 9rr 



>//o'- (3-4) 




But this implies that, just as in the DGP example,^ 

1/3 

(3.5) 

and not Tc = Hq^, is the distance smaller than which one expects deviations from Einstein 
gravity to be small, but larger than which significant deviations from Einstein must occur 
in the force law in order to reproduce the desired cosmic history as a modification of the 
Friedmann equation. Indeed, we may write our modified Friedmann equation in the very 
general form: 

H'^H^gix), (3.6) 
with the dimensionless quantity x defined as 

g{x) is such that g{x) —>■ x when x ^ 1, but g{x) substantially deviates from x otherwise, e.g. 
for the cosmological constant case g{x) — x+ft\, with x — Jl^(l + 2;)^ from Eq. (3.7) and 
and JIa the energy densities of matter and cosmological constant in terms of critical density 
at redshift z — 0. The modified Schwarzschild metric may then be read from Eqs. (2.13) 
and (2.14): 

E-'goo = g-' = 1 - r'H^g (^) ■ (3-8) 

where now x = rl/r^, using the definition of r* and Eq. (3.3). 



^In DGP gravity, is the distance at which scalar would-be radion modes become free to propa- 
gate, adding a Brans-Dicke type scalar to the existing gravitational interactions, where the strength 
of the Brans-Dicke coupling, u, depends sensitively on the background cosmology. 
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IV. GROWTH OF PERTURBATIONS 



Wc next wish to study the influence of the modified gravitational force law on the evo- 
lution of perturbations in the universe. Define the top-hat overdensity S{t) of a spherical 
mass of dust with mass M and radius r by 

M 
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(4.1) 



where p{t) is the background matter density (dust component only, i.e., not including the 
energy density of the dark energy). The conventional evolution of S{t) is governed by the 
following equation [21] 

4 



S + 2HS - 



31 + 5 



5^ = AnGp{l + S)S 



(4.2) 



The quantity H{t) corresponds to the background evolution. We wish to compare this 
relationship with that for the same scale factor evolution (and correspondingly, the same 
p{t) evolution), but with modified gravity. 

By exploiting the Birkhoff 's law constraint, one may compute the evolution of overdensi- 
ties by merely following the geodesies of spherical masses, without regard to physics outside 
the spherical mass itself. Note that this is not possible unless spherically symmetric config- 
urations respect the metric Eqs. (2.13) and (2.14). (DGP gravity in particular does not fall 
under this category of Schwarzschild-like metrics.) Using the geodesic equation as expressed 
by differentiating Eq. (2.14) with respect to t, we get 



1 d 



^xg'{x) 



(4.3) 



Note that this resembles Newton's second law, but is fully relativistic. Using some algebra 
and rewriting Eq. (4.1) as follows 



1 + 5 



X 



3Hl 
SnGp 



(4.4) 



where again x — one may use Eq. (4.3) to derive a new governing equation for 5{t) 



5 + 2H5 - 



4 1 



31 + 5 



\x{l + 5)g'{x{l + 5))-g{x{l + 5)) 



-3(1 + ml 



-xg [X) - g{x) 



(4.5) 



= The quantity a{t) is the background scale factor. Note 



where we define x 

that Eq. (4.5) is dependent only on the background evolution and makes no reference to 
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the mass M or the radius r*. Evolution begins with x ^ 1 and x decreasing with time. 
Deviations from the usual Einstein evolution occur when x ~ 1. 

Equation (4.5) is equivalent to Eq. (50) in [22], where it is derived assuming the continuity 
equation and the Friedman equation for fluctuations. It turns out that this is equivalent to 
assuming the validity of BirkhofF's theorem. As wc mentioned above, however, this approach 
is not justified when dealing with DGP gravity; basically, one cannot infer the evolution of 
a spherical perturbation from the evolution of the scale factor. 

In linear perturbation theory, we may simplify Eq. (4.5): 

5 + 2H5 ^A7rGp6[g'{x) + 3xg"{x)] . (4.6) 

li g{x) — X + Ai + Ai x^/^, where A\ and A2 are constants, we recover the usual scenario 
Eq. (4.2) and modified gravity yields an identical answer to the usual dark energy scenario. 
Choosing the time variable to be x, we see that 



3 

g{x) + -xg\x) 



% = \ b'(^) + 3:r^"(:r)] b . (4.7) 



Compare this expression with the corresponding one for evolution in a dark energy back- 
ground from Eq. (4.2) 



3 

9{x) + ^xg\x) 



^ = ^(5 . (dark energy) (4.8) 



It is interesting to note that, unlike Eq. (4.8), the case of modified gravity Eq. (4.7) has a 
decaying mode solution D^'^' (x H oc g^^"^ for arbitrary expansion histories. As a result of 
this, one finds that the growing mode D^'^ obeys 

z^r(.)4«;„'v'^/;jj(;^ (4.9) 

where D^'^ is normalized so that it scales like the scale factor at early times, x^ 1. When 
the expansion history can be described by a linear combination of dust, vacuum energy 
and curvature, g{x) — x -\- Ai-\- A2 x^/^, dark energy and modified gravity give rise to the 
same linear growth of density perturbations and Eq. (4.9) becomes the standard quadrature 
representation [23]. For a general expansion history, one cannot use Eq. (4.9) to find the 
growing mode in a dark energy background and must solve Eq. (4.8) instead. 

Let us now consider how the linear growth of perturbations can distinguish between 
dark energy and modified gravity models with the same expansion history. This is relevant 
because in the near future, planned experiments will measure the expansion history from 
2; = to 2; ~ 2 with very good accuracy (see e.g. [24,25]). For other tests of gravity using 
large-scale structure see e.g. [26,27]. To illustrate our results, we assume a modified Friedman 
equation given by 
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FIG. 1. Growth of perturbations as a function of redshift in dark energy (DE) models and 
modified gravity (MG) models with identical expansion histories, given by Eq. (4.10). The dashed 
line shows the ratio of density perturbation growth factors D^'^/DY^, the solid line shows the 
ratio of the velocity perturbation growth factors / f^^- 



where c is a constant, a = (1 — r2m)/2. This corresponds to "extra-dimensions-inspired" 
modified gravity in [20] when their a = 1 (which corresponds to an effective equation of 
state with ~ —0.7 for z between and 2). In all calculations we set Vt^n = 1/3 for 
simplicity. 

Figure 1 shows the results for the ratio D^^/ (dashed line) as a function of redshift 
z. The modified gravity perturbations growth slower, and the decline of this ratio as redshift 
approaches zero is a general feature, not restricted to Eq. (4.10), that can be understood in 
general terms. From Eqs. (4.6-4.8) we see that for the same expansion history, the difference 
in growth can be thought as coming from an effective gravitational constant G[g'{x) + 
3xg"{x)]. At high redshift g{x) ~ x, but at low redshift the constraint that the universe 
accelerates (a > 0) implies dlngf/dlnx < 2/3. If we model g{x) as a local power-law, 
g{x) ~ x", where n changes very slowly from n = 1 aX high redshift to < 2/3 when z ~ 1, 
we see that g'{x) + 3xg"{x) ~ n{3n — 2)g{x)/x, so as n becomes closer to 2/3, g'{x) +3xg"{x) 
becomes small. Although in practice one must take into account the time dependence of 
n, this illustrates why the growth of structure is slower in the case of modified gravity. As 
long as we are only concerned with linear perturbations around a cosmological background, 
this effective Newton's constant is wavelength independent. However, when perturbations 



c -|- vcF+x 



2 



(4.10) 



9 



become strong, and self-gravitation dominates cosmology, we will see that the effective G 
does acquire distance dependence (see Sec. VI). 

Figure 1 also shows the ratio / f^^, where / = dln£)+/dlna governs the growth of 
velocity fluctuations in linear perturbation theory, that is, the velocity divergence evolves 
as ^ = —HafS from the hnearized continuity equation (sec e.g., Ref. [21]). We see that 
slower growth leads to smaller time derivative, thus j^G^yDE ^-^so decreases as 2; — > 0. 
These deviations are detectable with precision measures of large-scale structure; for example, 
D+ and / can be derived from joint measurements of the redshift-space power spectrum 
anisotropy and bispectrum. The Sloan Digital Sky Survey should be able to probe these 
quantities with statistical errors of order a few percent [28]. 

So far we have discussed the hnear growth of perturbations. Equation (4.5) can be recast 
in a more useful way to study the non-linear evolution in perturbation theory. 



g{x) + -xg [x] 



dx 



di~ l + 5^d^^ ~ ^ 2x ^ ^1 n\ 

"(3n-2)y(") + 3xy("+^)l (4.11) 



A couple of points are worth stressing here. First, this equation becomes identical to that in 
the standard dark energy scenario, Eq. (4.2), only if g{x) — Ai + x, i.e. a curvature term is 
no longer degenerate as in the linear case. Expanded to second order, this equation can be 
used to compute the skewness of the density field (see also [22]); we have checked that for the 
same expansion history given by Eq. (4.10) the skewness changes by less than one percent 
in modified gravity compared to dark energy models. Only a rather large third derivative 
of g can induce an appreciable change in skewness from that in the dark energy case. Such 
extreme models have been studied in [22] (e.g. using Eq. (6.5) below with n — and q — 5); 
we have verified that in this case the skewness can change by up to 10% compared to dark 
energy models, and also D^^/D^ can be larger HiOXL unity (up to 1.2 at 2; = 0) due to a 
fast variation of g with x. 

Ultimately, unless the source of acceleration is a vacuum energy component, determina- 
tion of the expansion history of the universe plus the growth of structure should allow us 
to identify, through the divergence between Eq. (4.2) and Eq. (4.5), whether today's cosmic 
acceleration may be attributed to modified gravity or dark energy. 



V. EVOLUTION OF GRAVITATIONAL POTENTIALS: THE ISW EFFECT 

The microwave background provides another window to differentiate between dark energy 
and modified gravity through the late-time integrated Sachs- Wolfe (ISW) effect due to the 
decay in gravitational potentials at late times (see e.g. [29]). In order to assess how a 
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modified gravitational force law affects the ISW effect, we must find the evolution of the 
cosmological gravitational potentials, ^{t,X) and $(t,A), which are defined using the fine 
element 

ds"^ = [1 + 2'^ {t,X)]dt^ -a\t)[l + X)][dX'^ + X'^dn\ , (5.1) 

where a{t) is the background scale factor evolution. We are only interested in potentials and 
overdensities that are small and, therefore, only interested in linear perturbations around 
the cosmological background. 

We have the equations for the evolution of a hnear top-hat overdensity 5{t) of radial 
extent A*, Eq. (4.7), and we know the complete metric for that matter configuration: 

dX^ 



ds' = dt' - d\t) 



1-KA2 



+ x^dn 



(5.2) 



where = {1 + 6)a^, and k is the curvature parameter associated with the magnitude of 
the overdensity. Using r = aA* and integrating Eq. (4.3), we arrive at 

f^^il-KXl)-g;; , (5.3) 

where the constant of integration is determined by smoothly connecting a background cos- 
mological expansion outside A* with an overdense space inside A*. Then, one may deduce 
that 

2 

K, = Hlxg' {x)5 + -ddS . (5.4) 

One can confirm using Eq. (4.9) that k, is indeed a constant. 

One needs only identify a coordinate transformation taking the metric Eq. (5.2) into the 
form Eq. (5.1) and read off ^ and After some algebra we find 

^{t, A) = (l - ^] [inGpSig' + 3xg")] (5.5) 

GpSg'] , (5.6) 

when A < A* and the potentials vanish outside A*. By superposing top-hat linear overden- 
sities, we surmise that for general r-dependent 6 

V^* = AnGpa'^ [g'(x) + 3xg"{x)] 5 (5.7) 
= -A7rGpa^g'{x)S , (5.8) 

where the Laplacian is with respect to comoving coordinates. These expressions for the 
linear gravitational potentials are the generalization of the usual result for Einstein gravity. 
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FIG. 2. The angular cross-correlation function between galaxies and CMB temperature 
anisotropies. The solid line shows the dark energy case, dashed line corresponds to modified 
gravity with the same expansion history. 

around a specific cosmological background. One arrives at time-dependent effective Newton's 
constants (in general, different for \1/ versus $). At early times, Gefr G but deviates from 
the true Newton's constant significantly as H —>■ Hq. Note that this large discrepancy 
only applies to self-gravitation of linear density perturbations. For example, solar system 
manifestations of Newton's constant, in this context, are not linear and will only have small 
deviations from the true Newton's constant, as will be seen in the next section. 

We may now take Eqs. (5.8-5.7) and apply them to ascertain the ISW effect on the 
cosmic microwave background, which is proportional to the integral of $ — \1/ along the line 
of sight (see e.g. [29]). From Eqs. (5.7-5.8) we find, after a Fourier transformation, 



where A = —k'^/{d)TiGHpa^5Q) with k the comoving wavenumber and 5q the amplitude of 
density perturbations at some early time. The second term in Eq. (5.10), representing the 
time derivative of the effective gravitational constant, leads to an additional decay of the 
potentials that can be quite significant at low redshifts. 

We compute the cross-correlation between temperature anisotropies due to the ISW 
effect and galaxy fluctuations, which has recently been detected using the WMAP and 
SDSS datasets [30,31]. The angular cross-correlation function WTgiO) can be written as 




(5.9) 
(5.10) 



MG 



(1 - /^^)(^' + Ixg") + \{^xg" + ?,x'g"')\ D 



iMG 
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TT^ r r die 

w^fiO) = 3T,n^bi2n f^ J dz^ - /) w,iz) J —P{k)Mk9x) , (5.11) 
for the dark energy case, whereas for modified gravity 

w^^'^ie) = 3Ton^b{27ry^ J dz^ 0% [(1 - f){g' + ^-xg") + ^x{bg" + 2>xg"')\ w,{z) 

X j ^P{k)Mk9x) , (5.12) 

where we have assumed the small-angle approximation, and Tq is the CMB temperature, b 
is the hnear bias of LRG galaxies, Wg{z) is the galaxy selection function, P{k) is the dark 
matter power spectrum at z — 0, Jq is a Bessel function, and x(z) is the comoving distance 
as a function of redshift z. We assume, following the results in [31] (see their Fig. 3), that^ 
b — 5.47, and use their selection function for the z — Q.AQ sample. 

Figure 2 shows the results, again using the expansion history given by Eq. (4.10). We see 
that the MG angular correlation function is suppressed by a factor of about two compared 
to the case of DE, for the same expansion history. This anomaly is a reflection of the 
order-unity anomalous effective Newton's constant seen in Eqs. (5.8-5.7) when the Hubble 
parameter is near today's value. Although the ISW anomaly is not yet detectable, it should 
be so in the near future with the completion of the SDSS survey or other future probes of 
structure at intermediate redshifts. 

The suppression of the late-time ISW effect that results from MG can help explain, at 
least partially, the low amplitude of the CMB power spectrum at low multipoles as confirmed 
by WMAP [32]. It is worth emphasizing that these results depend on the specific prescription 
of exploiting Birkhoff's law for understanding modified gravity from cosmological physics, 
and it will be important to check whether the large ISW discrepancy carries over to more 
general prescriptions of modified gravity. 

VI. ORBIT PRECESSION 

Wc have seen that modification of the Friedmann equation leads to a modification of the 
gravitational force law. This modification of the force law may, in turn, lead to small-but- 
detectable corrections at distances much smaller than today's Hubble scale. Indeed, under 
some circumstances (e.g., if no physics other than that deduced from cosmological evolution 
emerges at small scales) these alterations of cosmology on the largest observable distances 



^This choice is a matter of convention here since it scales both predictions by the same amount. In 
practice, by measuring the angular bispectrum of LRG galaxies one can determine h, which due to 
the results in the previous section should be almost independent of whether DE or MG is present. 
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may cause observable deviations from known physics at even solar system scales. Let us 
elaborate. 

The precession A0 of the perihelion per orbit in the background of a metric of the form 
Eq. (2.2) may be determined in the usual way: 



= y dr^ 



E^-gooil + J/r^) ' 



(6.1) 



where E = N'^T, J = r'^cj), and dot refers to differentiation with respect to proper time. For 
a nearly circular orbit with a metric of the form Eq. (3.8), one may compute the precession 
rate: 



^ -f^O 1/2 



9 - l^g' 



(6.2) 



with again, x = The leading contribution to orbit precession from the altered metric 
Eq. (3.8) comes from the simple alteration of the Newtonian potential. 

We are particularly interested in orbits whose radii around a central body of 
Schwarzschild radius Tg are much smaller than = {r-g/ HqY^^. Then, corrections from 
the modified gravity are small. One may represent the function g{x) as 

g{x) =x[l + 6g{x)] , 

with Sg{x) < 0(1). Then, Eq. (6.2) reduces to 

{Sg' + f 5^" 



dT ^ V2 



(6.3) 



(6.4) 



to leading order in 6g. Recall that here x ^ 1. 

Take as an instructive example the form of the modified Friedmann equation found in 
Cardassian models [19]. This can be written as 

1/9 



g{x) ^x[l + cr-«('-") 
where n and q are parameters of the modification and where 

l^l + (l + Zeq)=^^a-"), 



1-n 



(6.5) 



(6.6) 



The quantity z^q^ is the redshift at which the two terms inside bracket in Eq. (6.5) are equal, 
a quantity of order unity. Then, using Eq. (6.4), the anomalous orbit precession rate is 



dT ^ 



^(n-l)[3.(n-l).l] y 



(6.7) 
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When q{l — n) — |, this expression is independent of the radius of the orbit, or the mass 



of the central body, and the anomalous precession rate is proportional to today's Hubble 
parameter. Ho ~ 10 //as/year. Such a precession rate is on the threshold of detection by 
precision ephemeris measurements of the inner solar system, particularly with intriguing 
developments this decade coming from two Mercury-bound missions (BcpiColombo and 
MESSENGER) as well as improvements in lunar ranging observations [14,15,33-36]. When 
q{l — n) 7^ |, there is a relevant distance-dependence for the anomalous precession rate, 
where the governing length scale is once again r*. Since we are primarily concerned with 
orbits such that r -C r*, the dimensionless distance factor in the anomalous precession rate, 
j'j,3^j,3jq(i-n)-i^ either be huge or tiny. Thus, with solar system constraints in mind, the 
parametric range where q{l — n) ^ ^ can already be ruled out. However, for g(l — n) ^ 
no solar system test are likely to discover discrepancies based on anomalous orbit precession 
in the foreseeable future. 



In this paper we showed how modifying gravity to effect the observed late-time cosmolog- 
ical acceleration at scales of today's Hubble radius, H^^, can lead naturally to corresponding 
modifications of gravitational interactions at scales much shorter than H^^. Indeed, by pre- 
suming that the new gravitational physics obeys a limited version of Birkhoff 's law, we were 
able to derive the precise form of the modifications to Newton's law of gravitation at short 
(sub-cosmological) distances. We then showed that an observer in the gravitational field of 
a central source whose Schwarzschild radius is Vg, experiences substantial deviations from 
the usual Schwarzschild metric at all distances greater than approximately 



For many models these deviations will be measurable through observation of orbital pre- 
cession of solar system objects in the coming decade. We also discussed the evolution of 
density perturbations and showed that, unless the acceleration of the universe is driven by 
an effective vacuum energy, simultaneous measurement of the expansion history and growth 
of large-scale structure can be used to distinguish modified gravity from dark energy. In 
addition, the cross-correlation of galaxy distribution and the cosmic microwave background 
temperature anisotropy can detect anomahes in the late-time integrated Sachs-Wolfe effect 
caused by modified gravity. Such measurements will be available imminently. 

It is instructive that these results arc identical to those found for the braneworld theory 
of Dvah, Gabadadze, and Porrati (DGP), even though DGP gravity does not respect any 
dynamical version of Birkhoff's law. The correspondence between the scales of departure 



VII. CONCLUDING REMARKS 




(7.1) 
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from Einstein gravity in DGP and the Birkhoff 's Law theories seems not to be a coincidence. 
Indeed, one suspects that it is quite general. Imagine cosmology at extremely late times, 
when all matter surrounding a particular gravitational source is swept away. Then, if one 
believes that an isolated, central source has a (quasi)static metric description, it should be 
Schwarzschild at short distances and deviate from Schwarzschild at large-enough distances. 
How large? Since this metric must still encode the cosmology within it, i.e., test observers at 
large distances from the source should recede from the source in the manner dictated by the 
given late-time (accelerating) cosmology, the empty space metric must include a repulsive 
force at distances where cosmological flow overcomes the local gravitational binding of the 
central source. Thus, we expect substantial deviations from Schwarzschild (in the form of a 
repulsive force) at distances C(r*), Thus, it is not unreasonable to expect general modifica- 
tions of cosmology at today's Hubble scale to inevitably affect local gravitational interactions 
at distances governed by r*; it is only the precise functional form of the deviations that will 
vary from model to model. 
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